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Abstract 

The dynamics of linear perturbations in a radiatively cooling dusty plasma 
is considered, with the charge of both dust (Zd) and plasma (Z p ) components 
being allowed to vary with their densities. It is shown that in the long- 
wavelength limit corresponding to the characteristic cooling length, when the 
plasma can be treated as quasineutral, the presence of dust particles changes 
the criteria for radiation instability, regardless the charging process of the dust 
particles. In particular, the condensation (isobaric) mode is shown to be stabi- 
lized (destabilized) if in the equilibrium, the relation between densities of the 
dust nd and plasma n under the quasineutrality condition, [d In n^/d In n) q < 1 
(> 1) is satisfied, while the isentropic mode is stabilized (destabilized) when 
the opposite inequalities take place; the isochoric mode is unaltered. Numeri- 
cal estimates show that these effects can be important in hot phases (T ~ 10 6 
K) of the interstellar plasma, and in tokamak plasma near the walls. 
PACS numbers: 95.30.Qd, 95.30.Lz, 52.35.-g, 52.35Py 
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I. INTRODUCTION 



It is well known (see, for example, Refs. 1, 2, 3) that the presence of charged dust particles 
in plasma changes its dispersion properties qualitatively, and that the dust itself can give 
rise to new plasma instabilities. From intuitive point of view it seems possible that dust 
must be important in dynamics of radiatively cooling plasmas, not only because it absorbs 
radiation, but also due to its ability to transform thermal energy of the plasma into radiation 
in inelastic collisions with electrons and ions. Apparently, it can be of great importance in 
high temperature (T >, 10 6 K) cosmic and laboratory (in particular, close to tokamak walls) 
plasmas, where frequent collisions with the ions maintain dust particles sufficiently hot, so 
that they re-emit thermal energy efficiently 4 . The present work is aimed to analyze the effects 
of charged dust particles on the conditions for the radiation instabilities. It is important to 
note that this issue has been addressed recently in 5 ' 6 , however in these papers dust particles 
were considered only from the point of view of electrostatic interactions, while the cooling 
processes were treated as independent on the presence of dust. In our study we explicitly 
assume that dust particles contribute substantially to the net cooling rate, which brings to 
the system qualitatively new dynamical features. For the sake of simplicity we will consider 
here only highly ionized plasmas, when the effects from neutral particles can be neglected. 
In Sec. II. A we briefly discuss how important can be dust as a cooler of the plasma. In 
Sec. II. B we formulate a simplified two-fluid model of a collisional and radiatively cooling 
dusty plasma. In Sec. III. A we describe the steady state equilibrium, with the relevant 
parameters estimated in Sec. III.B, while in Sec. III.C linearized equations are written. 
In Sec. III.D the dispersion equation and its solutions are given and criteria of thermal 
instability in different wavelength ranges are discussed. Sec. IV summarizes the results. 
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II. BASIC EQUATIONS 

A. Radiatively cooling dusty plasma 

In an optically thin plasma dust particles, if present, besides their dynamical effects through 
electrostatic interaction with plasma particles, can also work as a cooling agent via the 
collisional transfer of thermal energy of plasma component into infrared radiation. If one 
assumes, as an idealization, that dust particles emit as black body, the energy loss rate due to 
dust particles varies as oc T d , where T d is temperature of the dust. For realistic particles this 
dependence can be steeper. Detailed calculations of the relevant cooling rate are made in Ref. 
4. They show that for physical conditions and abundance of dust typical for the interstellar 
medium (ISM) , radiation losses from dust particles dominates cooling from the line emission 
of heavy elements by two orders of magnitude in the interval T = 10 6 — 10 10 K 4 . One should 
stress, however, that at temperatures T ~ 3 x (10 7 — 10 8 ) K dust particles are efficiently 
destroyed in collisions with the ions and electrons 7 , and therefore only in temperature range 
between 10 6 K and 3 x 10 7 K dust contributes sufficiently to the net cooling. Apparently, 
dust cooling can also be important in tokamak plasma near the walls, where temperature 
can lie in this interval, or in regions of MARFEs (multifaceted asymmetric radiation from 
the edge) formation if radiating impurities are introduced into the plasma by sputtering. 

B. Two fluids 

We will consider a simplified two-fluid approach, with the dust charge instantaneously 
(on times much shorter than other characteristic time scales of the system) settled on the 
equilibrium Z d = Z d (n,n d ,T), where n and n d are the plasma and dust densities, T is the 
plasma temperature. To outline the main physical effects, we will study the simple case of 
potential perturbations in the planar geometry. The full set of equations in this case is 

d t n + d x (nv) = 0, (1) 
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pd t v + pvd x v = -d x p - eZ p nd x (j) - f Q (v — u), 



(2) 




£(T, n, n d ) 



(3) 



d t n d + d x {n d u) = 0, 



(4) 



p d d t u + p d ud x u = -eZ d n d d x 4> + fo(v - u) 



(5) 



p = k B nT, 



(6) 




47re(Z d n d + Z p n), 



(7) 



where v, u are velocities of the plasma and dust particles, respectively, fo [g cm~ 3 s _1 ] is the 
coefficient describing the friction between dust and plasma, p = mn and p d = m d n d are the 
mass densities of the plasma and dust components, p, the gas pressure, 0, the electrostatic 
potential, Z d , the dust charge (negative for negatively charged dust), Z p , the mean plasma 
charge defined as Z p n = (ZiU — n e ) = {Z^ — x)n with Zi being the ion charge, and x = n e /n; 
note that for positively charged dust particles x > Zf, here an adiabatic index 7 = 5/3 has 
been assumed and the net cooling function £ explicitly depending on number density of dust 
particles; the dust component is assumed kinetically cold, so that the dust pressure term in 
Eq. | is omitted; d a and d a denote partial and Lagrangian derivatives with respect to a. In 
general, dust particles contribute both to heating {e.g. via photoemission of electrons from 
grain surface, see for recent discussion Ref. 8), and cooling (via collisional transfer of kinetic 
energy of plasma particles into heat of dust grains and subsequent re-emission in, as a rule, 
infrared range 4 ). It is obvious that the charge separation are important only on the spatial 
scales comparable (or shorter) to Debye length, on longer lengths plasma can be treated as 
quasineutral (see, below). 
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III. DISPERSION RELATION 



A. Unperturbed state 

We assume that in an unperturbed state radiative cooling is balanced by energy input 
from external sources (e.g. hard emission, Ohmic heating etc.), so that 

C = 0, (8) 

[note that in most cases the net cooling function can be presented in the form C = 
L(T,n,rid)n — T(T,n,nd)n, where L(T,n,nd)n is the radiative cooling rate, T(T, n, rid)n, 
the heating rate due to external energy sources]. We assume also that unperturbed plasma 
is quasineutral 

Zdondo + Z p0 n = 0. (9) 

We will further assume that the dust and plasma charges adjust the corresponding equi- 
libria instantaneously Z d = Z d (n,rid,T) and Z p = Z p (n,n d ,T), i.e. that the charging 
characteristic times are shorter than other relevant times. 



B. Evaluation of the relevant numbers 

The problem is characterized by three time scales: thermal, frictional, and electrostatic. 
Thermal time scale is written as 

70 = —j—, (10) 

here L is the radiative cooling rate taken at the equilibrium, Boltzmann constant. 
Frictional times for the plasma and dust particles 

Jo Jo 

respectively, and electrostatic time (the inverse plasma frequency) 



All characteristic times (frictional and electrostatic) relevant to the dust component are 
connected with the corresponding plasma times by the intrinsic constants 5 = n d o/n , the 
concentration of dust particles, and p = m d n d0 /mn = pdo/po, the mass density ratio of the 
dust and plasma components: 

T f4 = fJ>T f , r e4 = «J-V -1/2 r e . (13) 

It is readily seen that two dimensionless parameters appear in the problem 



r k B T /47re 2 n r k B T f 

k= — = — — \ , andz/= — = — , 14 

r e L V m r f L p 

which characterize the role of electrostatic and frictional forces in dynamics of radiation- 
condensation instability. 

C. Linearized equations 

The full set of linearized dynamical equations for a nonmagnetized dusty plasma is 
written as 

d t n + n d x v = 0, (15) 

mn d t v + d x p + eZ p0 n d x (p + f (v - u) = 0, (16) 

\d t p - l^d t n + (d n C)n + (d nd C)n d + (d T C)T = 0, (17) 

d t n d + n do d x u = 0, (18) 

m d n do d t u + eZ do n do d x <i) - f (v - u) = 0, (19) 

p n T , , 

Po n T 
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dlcf) = -47ie[Z d0 n d + (d n Z d )n do n + (d T Z d )n dQ T + (d nd Z d )n do n d + 

Z p0 n + (d n Z p )n n + (d T Z p )n T + (d nd Z p )n n d ], (21) 

where the unperturbed variables are given with subscript zero, a , while the perturbations 
are without subscript, a. 

In order to evaluate the interrelation between electrostatic and radiative effects it is 



useful to write the linearized Poisson equation 21 in a non-dimensional form 



«9|0 = -K 2 [5{Z dG n d + {d n Z d )n + (d nd Z d )n d + (d T Z d )T} + 

{d n Z p )n + {d rid Z p )n d + {d T Z p )f + Z p0 n], (22) 

here all variables with bar are normalized to unperturbed values: a = a/a^ t = t/r Q , 
£ = x/t c , where c is the isothermal sound speed, the potential is normalized to mcp/e, n 
is defined above. 

For typical interstellar or tokamak hot plasma at T ~ 10 6 K, k can be estimated as 

2 3 ■ 10 30 

k ~ , (23) 

n 

It is readily seen from these estimates that motions on electrostatic and radiation scales can 
be separated due to a huge difference in corresponding times and lengths. This simply means 
that electrostatic forces acting between the charged dust and plasma components keep dusty 
plasma quasineutral on much shorter time scales than the relevant hydrodynamical scales. 
Formally it can be expressed in the form of perturbation procedure applied to the above 
dimensionless linearized equations with dynamical variables as a series 

a = a + kTV + .., (24) 

and potential as 



K 



+ + (25) 



Thus in what follows we turn to the dimensional variables and equations, with the quasineu- 
trality condition instead of solving Poisson equation. 



D. Dispersion equation 



In this framework the quasineutrality equation can be written in terms of n d 

n d = - [Z d o + {d nd Z d )n d Q + (9 nd Z p )n ] _1 [(9n-^d)ndo + Z p0 + (d n Z p )n ]n + 

-[Z d0 + (d nd Z d )n dQ + (d^Z^no^KdrZ^rido + (d T Z p )n ]T. (26) 

From Eq. (]2§) it follows that {d T n d ) q = -[Z d0 + (d nd Z d )n d0 + (d nd Z p )no\- 1 [(d T Z d )n m + 
(dTZ p )n Q ], which will be taken as zero in order that the quasineutrality to hold for purely 
isochoric perturbations; here the subscript q denotes that the derivative is taken over the 
quasineutrality state. 



Combining Equations of motion [16| and ^ in order to eliminate undefined we arrive 
finally for perturbations in the form a oc exp(ttt + ikx) to the eigenmatrix M for the 



eigenvector (n , T ,v ,u 

( 

M = 



n 



no 



ikriQ 



\ 



ik^ p tt p do tt 



■^tt + £ n + A£ nd lftt + C T 







V 



(27) 



Att ikn d0 

where A = -[Z d0 + (d nd Z d )n m + (d nd Z p )no\- l [(d n Z d )n do + Z p0 + (d n Z p )n ] = (dn d /dn) q . 

Finally the compatibility condition det(M) = leaves the dispersion equation 
2 T'o „ , , o 5,oo„ 2 



fi 3 + \^C T tt* + ^cltt + \k\\ 
6 po o 6 



—Ct (£ n + AC n ) 

Po Po 



0. 



(28) 



where c d = co/a/1 + /x, c being the isothermal sound speed in a plasma gas. It is seen that 
this equation is identical to the standard dispersion equation 9 , when all terms correspondent 
to dust component are omitted. 



Dl. Isobaric mode 

In general, there are three solutions of Eq. |28|. One corresponds to the so-called isobaric 
(condensation) mode, which formally can be obtained from Eq. putting k — > oo and 
Im(tt) = (see discussion in Refs. 9, 10), i.e. 
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The condition for instability of this mode is also straightforwardly follows from the Hurwitz 
criterion that the last coefficient of the third order polynomial in Eq. ^ is negative, i.e. 
n (C n + AC nd ) - T C T > 0. 

For the cooling function in the form C = A(T)n 2 + Ad(T)nrid — Tn, one can write the 
growth rate as 



5p 



1 + 



' d In rid 
d\nn 



- 1 

g 



Vd 



d In C c 
dlnT 



(30) 



where C c = A(T)n 2 + Ad(T)nnd, r/d = Ad(T)nn,d/ £ c is the relative contribution of dust 
cooling. It is readily seen that for strongly dynamically coupled ions and dust when rid = 
5n, electrostatic forces do not affect radiation instability. When ( d In rid/d In n) q < 1, i.e. 
quasineutrality and dust charging requires that the dust particles escape compressed regions, 
the condensation mode is stabilized because the corresponding condition 

din rid 



d\nC c 
~dTn~T < 1 



dlnn 



1 



Vd, (31) 



fulfills in a narrower temperature interval than when electrostatic effects disappear. Con- 
trary, if quasineutrality and dust charging requires an enhancement of dust particles in 
compressed regions, i.e. (dlnn^/ d\nn) q > 1, the condensation instability is enhanced. Qual- 
itatively it can be understood as due to an enhancement of the net cooling rate associated 
with dust cooling. 

D2. Isentropic (acoustic) mode 

The other two modes in the short- wavelength limit correspond to overstable acoustic motions 
with the wave numbers ±k and the growth rate much smaller than the wave frequency: 
|i2e(fi)| |/m(f2)| 9,10 . In order to find the solution in this case we will expand it as 

/5X 1 / 2 

Q = ±i - kc d + Or, (32) 
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where Q r is the growth rate of the amplitude of acoustic waves, for which we obtain to the 
first order to k 



15p 



TqC t + ^n (C n + AC nd ) 



(33) 



or 



2£ c (d\n£ c 3 
T5p~o{ d\nT + 2 



1 + 



d In rid " 
dhxn 



- 1 



Vd 



(34) 



It is seen from fl3~4] ) that contrary to the isobaric mode, isentropic perturbations are destabi- 
lized if (d In rid /d In n) q < 1, since the restrictions on the cooling rate C c is here weaker than 
in the absence of electrostatic forces. Instead, isentropic perturbations are stabilized when 
(d In rid/d In n) q > 1. Such a behavior is qualitatively clear because the instability of isen- 
tropic mode is physically connected with overheating of adiabatically compressed regions, 
and as soon as dust cooling decreases in compressed regions when (d\nrid/d\nn) q < 1, it 
results in an additional overheating. 



D3. Isochoric mode 

The instability criterion for the isochoric mode corresponding to the long-wavelength limit 
k — ► 0: Ct < 0, remains obviously unaltered by the electrostatic forces and additional dust 
cooling, because this limit describes perturbations with T/0 and rid — n — 0. 



V. SUMMARY 

In the present paper we considered radiation instability in an optically thin dusty plasma 
with dust particles contributing substantially to the radiative energy losses. We have shown 
that combination of the two effects: electrostatic interaction of the dust and plasma com- 
ponents and their tight dynamical coupling (expressed in quasineutrality) , and additional 
energy losses connected with dust particles, strongly modifies the criteria for the instability. 
This influence differs for different instability regimes: the isobaric (condensation) mode is 
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shown to be destabilized in dusty plasmas if the dust particles are excessively attracted in 
the regions of gas compression, (d\nnd/d\nn) q > 1, and stabilized in the opposite case; in 
the short-wavelength limit corresponding to the isentropic (adiabatic) mode effects from the 
presence of dust work differently: they quench the instability when (d\nnd/d\nn) q > 1, and 
destabilize it in the opposite case; the condition for the instability remains unaltered in a 
long-wavelength range where isochoric mode dominates. 

In a simplest case of a dilute hot plasma with dust particles charged predominantly 
collisionally, i.e. j p — > 0, (d\nrid/d\nn) q approaches unity from below as (see Appendix) 

\ll^)= l -\mj (35) 

and therefore the isobaric mode is quenched, while the isentropic is weakly enhanced. 

As j p approaches the critical value j pc = (m i /m e ) 1 '' 2 — 1, when dust particles change 
the charge from negative to positive, (d\nrid/d\nn) q — > — oo at j p = j pc — 0, and 
(d\nnd/d\nn) q — > +oo at j p = j pc + 0, (see Appendix). Therefore, from the left side of 
the critical point, j p = j pc — where the charging is dominated by collisions, the isobaric 
mode is strongly quenched, while the isentropic is strongly enhanced. Contrary, from the 
right side of the critical point, j p = j pc + where the grains are charged predominantly by 
photoionization, the isobaric mode is enhanced, while the isentropic one is suppressed. 

One can therefore conclude that in hot and dilute astrophysical plasmas embedded in 
a strong radiation field, where radiation dominates dust charging, the thermally unstable 
isobaric mode is enhanced. Under the conditions of an exposing radiation field with an 
intensity well below the critical value (such as in astrophysical plasmas far from strong 
radiation sources and dusty plasma sheaths near the tokamak walls) when the collisional 
processes are the principal charging mechanisms, isentropic mode is destabilized, while the 
isobaric one is quenched. 
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APPENDIX 

Generally, the equation for dust charge is written as (e.g. Refs. 1, 7, 11, 12) 

/ m \ 1 / 2 

3 P +(l-y)(l-Qy) = ( — ) (1 + Qy)exp(y), (Al) 

where y = e((p s — (4>))/kT, (fi s is the grain surface potential, (0) is the average potential 
in plasma, j p = I p /(eViCi 2 n), the photoionization rate for the grains in units of the ion 
charge flux on the surface, Q = AnX^naC, C = a(l + cl/Xd) is the capacity of a dust grain, 



A d = yJkT/[4:ir(n e + n)e 2 ], the Debye length; for a dilute dusty plasma, ATca 3 n d /3 1 and 
a X D , Z d = kTay/e 2 and C = a. It is readily seen from Eq. |AT] (see discussion in Refs. 
1, 7, 12) that at j p — > dust charge is negative, Z d < 0, however, when radiation flux 
increases (or n decreases) dust grains become charged positively 12 . The dependence n d (n) 
along the quasineutrality condition can be found as the solution of Eq. along with the 
quasineutrality condition ^| written here as 

l-z = --C, (A2) 

yo 

where for a dilute plasma y = e 2 / (k B Ta), and ( = n d /n. With these notations Q = 
£/[(!+ x)yo\. Equations [Al] and |Al| must be solved with a complementary equation for the 



balance of electrons 

' rru n x l 2 



j P C - — - (1 + Qy) exp(y) + b e x a x + j p x a , (A3) 



which we assume to be dominated by photoemission of electrons from dust grains, pho- 
toionization of neutrals by collisions with thermal electrons, photoionization of neutrals by 
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external radiation, radiative recombination, and sticking of electrons on dust grains. Here 
a r = a r /(4irvia 2 ), a r is radiative recombination coefficient, b e = /3 e /(47TVja 2 ), (3 e is the rate 
of collisional ionizations of neutrals, j p = l p j (vid 2 n) is the photoionization rate of neutrals, 
%a = n a /n, their fraction. 

Assuming that |1 — x\ <C 1, and taking into account that for hot plasma (T > 10 6 K) 



Or, b e ~ 
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-12 



10 11 <C Jmi/m e , one can reduce Eqs. |A1~ HA3| to the following system 



{1 ~ x)y _ (mj/m e ) 1/2 

jpC ~l~ jp-Ea 



c 



X 



(mi/m e ) l l 2 C 

which has for i a C 1 one physically meaningful root 



C ~ [(mi/m e ) 1/2 - j p ]y (m e /mi) 1/2 



In 



{mj/m e ) l/2 
1+J P 



(A4) 
(A5) 

(A6) 



The other root vanishes at j p and thus is unphysical. It is seen that at j p = j pc — 1 = 
[(m^/me) 1 / 2 — 1]— the solution £ — > +oo logarithmically, while from the right side j p = j pc +l 
as ( — > — oo. The derivative 



dn d 
dn 



n (|) +c ^- ( [ln w 1/2 i-^ +0 f [ln 



(mi/m e ) 1/2 / 
1+J P 



(A7) 



1+J P 

and hence always goes to — oo at j p ~ j pc . Note, that at this point x ~ 1 — (me/m^) 1 / 2 , so 
that |1 — x| 1 holds. The logarithmic derivative 



rz,i tin 



In 



{mi/m e ) 1 / 2 



+ 1, 



(A8) 



and goes to — oo from the left side of the point j p = (mi/me) 1 / 2 — 1, and to +oo from the 
right side. 

It is readily seen that in the limit j p — > 0, when charging is predominantly collisional, 

and a < Ad 

c^y /2 fe^^, (A9) 



and 



/ [^(mi/me) 1 / 2 - j p ] 

d\nn d \ ^ _ /77tj\ x / 2 
dlnn / \m P 



(A10) 
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